Unique Solutions to Hartree-Fock Equations 
for Closed Shell Atoms 



Marcel Griesemer and Fabian Hantsch 
^ , Universitat Stuttgart, Fachbereich Mathematik 

CN ! 70550 Stuttgart, Germany 

O 

^ , Abstract 

cn 

04 ' In this paper we study the problem of uniqueness of solutions to the Hartree 

and Hartree-Fock equations of atoms. We show, for example, that the Hartree- 
, Fock ground state of a closed shell atom is unique provided the atomic number Z 

is sufficiently large compared to the number N of electrons. More specifically, a 
<^ . two-electron atom with atomic number Z > S5 has a unique Hartree-Fock ground 

, state given by two orbitals with opposite spins and identical spatial wave functions. 

This statement is wrong for some Z > 1, which exhibits a phase segregation. 

^ '. 1 Introduction 

0^ 

' The Hartree-Fock method is widely used in quantum chemistry for approximate elec- 

l/^ ■ tronic structure computations [11]. In the simplest case the state of the electrons is 

fSj . described by a single Slater determinant and one seeks to minimize the energy by vari- 

ation of the one-electron orbitals. This is done by some self-consistent field algorithm 
such as the Roothaan, or the level-shift algorithm [5]. While there has been remarkable 
progress recently in the analysis of convergence properties of these algorithms [U S] , it 
^ I is still poorly understood what is actually being approximated. To some extent this is 

' due to our ignorance about the set of critical points of the Hartree-Fock functional: we 

^ know that the Hartree-Fock functional for a neutral atom has a minimizer and infinitely 

many other critical points [161 117j . but uniqueness of the minimizer, e.g., is not known 
even in cases where it is expected. Neither is it known whether distinct methods for 
finding critical points lead to distinct critical points. Our goal is to give answers to 
such questions. 

In this paper we establish existence and uniqueness of solutions to the Hartree-Fock 
equations for positively charged atoms with prescribed filled shells. We consider both 
restricted and unrestricted Hartree-Fock theory. In the unrestricted case a filled shell 
refers to a set of qn'^ electrons which means that we take the number of all electrons, 
N, of the form N = q X^^^^ where 1 < ni < n2 ■ ■ ■ < ris and where q denotes the 
number of spin states. In the restricted case our notion of a shell is the usual one 
and hence all values of N that occur in noble gas atoms are admissible. For atoms 
with partially filled shells our uniqueness results will not hold. Our method is based 
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on a perturbation argument which exploits the fact that in the hmit of large atomic 
numbers, Z — )• oo, the electron-electron interaction energy is negligible compared to 
the total energy. This forces us to choose Z much larger than N ^ but thanks to a 
novel technique for comparing spectral projections, values of Z as small as 35 can be 
handled in the case of two-electron atoms. Also, we provide an example which shows 
that Z > A'' — 1 is not sufficient for our results to hold in general. As a byproduct of our 
methods we obtain uniqueness of the minimizer of the A^-electron Hartree functional 
for sufficiently large Z. 

The Hartree-Fock equations for an atom with atomic number Z > and A^ electrons 
are the set of equations 



uT r^^(^) + ^2^ / u~Tj dy = eiipi{x) (1) 

12^1/ fc^i^ F y\ 

for A^ functions tpi, . . . ,ipN G -L^(M^ x {1, . . . , q}) subject to the constraints 



ipi{x)ipj{x)dx = 5ij (2) 

and real numbers ei,...,e7v- Here and henceforth x,y,... denote elements (x, s) of 
M.^ X {1, . . . ,q}, J dx denotes integration with respect to the product of Lebesgue and 
counting measure, and |x — y| = |x — y|. Of course g = 2 for (spin-1/2) electrons but 
for later convenience we allow arbitrary g G N. We have chosen atomic units where 
h = 1, the mass m of the electron equals 1/2 and the Rydberg energy equals 1. 

The following theorem is our main result on critical points of the unrestricted 
Hartree-Fock functional. An analog result on a restricted Hartree-Fock functional is 
given in Theorem 14.11 

Theorem 1.1. Let q, N, s and ni, . . . , be positive integers with 1 < rii < n2 < ■ ■ ■ < 
Ug and N = qYlk=i ^1- Suppose that 

Z > -^(20Ar + 8V2iV), A, = n^^ _ (n, + l)-2. (3) 
Then the Hartree-Fock equations ([1]) have a solution ipi, . . . ,ipN with 

ei, . . . , ETV G ULi [ - Z^/nl -Z^/iuk + 1)^), (4) 

and the space spanned by ipi, . . . ,lpn is uniquely determined by condition (jl]). The 
orbitals tpi may be chosen in the form 

yn;mo-(x,At) = •^"1^1^^ V ;n^(x)(5cr,^ (5) 

where n G {ni, . . . , Us}, 0</<n — 1, —l<m<l,a^ {1, . . . , q}, and each of these 
quadruples {n,l,m,a) occurs exactly once. 
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The Hartree-Fock equations ([T]) are equivalent to the Euler-Lagrange equations of 
the Hartree-Fock functional 



N 

£^^{ipi,...,ipN)='^ J \Vipk{x)\'^ - j^\^Pkix)\'^ dx 



k=l 

|2|,„. f,,\\2 



I - Yl [ ^"^'^^^^ \'fik{y)\ - ^i{x)ipi{y)^pk{x)ipk{y) ^^^^^ 



i<k 



where ipi,...,(pi^ G L^(M'^ x {1, . . . , g}) are subject to the constraints ([2]). In fact, the 
Euler-Lagrange equations of ([6]) are equations of the form (JT]) with ei^pi replaced by the 



more general term Yli^=i ^ij^ji coefficients Ay being Lagrange multipliers. These 
generalized Hartree-Fock equations, as well as the Hartree-Fock functional and the con- 
straints dl]) are invariant with respect to transformations {ipi, . . . , ip^) i— )• (^ipi^ . . . , (^tv) 
of the form ipi = Y^f=i Uij^pj where (Uij) denotes any unitary N x N matrix. By choos- 
ing this matrix to diagonalize the self-adjoint matrix (Xij), the equations ([1]) emerge. 
This also means that the Hartree-Fock functional only depends on the A^-dimensional 
subspace of L^(M^ x {1, . . . , g}) spanned by ...,(/? at or on the orthogonal projection 

N 

P = Y,\Vk){^k\ (7) 

k=l 

onto this space. As a function of P the Hartree-Fock functional is quadratic and its 
domain can be extended to a convex set without lowering the minimum [HIE]. Nev- 
ertheless the Hartree-Fock functional is not convex due to the presence of the negative 
exchange term - the second term in the numerator of ([6|). The convex functional 
obtained by dropping the exchange term is called reduced Hartree-Fock functional. 

We are now in position to derive Theorem 11.21 below, on the uniqueness of the 
minimizer of the Hartree-Fock functional ([6]). Suppose we choose ni = l,n2 = 2, . . . 
in Theorem 11.11 Then, by Proposition 12.21 for ^ sufficiently large the condition 
becomes equivalent to the statement that ei, . . . ,eN are the lowest eigenvalues of 
the Fock operator, which is the linear operator acting on Lfi on the left hand side of 
([!]). This condition on ei, . . . ,e]y is satisfied for any solution of ([T|) associated with 
a minimizer of the Hartree-Fock functional |16j . Hence the Theorem 11.11 implies the 
following theorem with the exception of the bound on Z. 

Theorem 1.2. Let s, q and N be positive integers with N = and suppose 

that Z > Aj^{12N + A^m - 4) where A, = s'^ - {s + 1)"^ Then the minimizer 
of the Hartree-Fock functional is unique in the sense of a unique projection ([7]). It 
may be represented by N orbitals ipi, . . . ,ip]\f of the form ([5]) satisfying the Hartree- 
Fock equations ([T]). In particular, the density ^^j^ is unique and spherically 
symmetric. 

For N = q = 2 and Z > 35 the Theorem 11.21 savs that the Hartree-Fock functional 
has a unique minimizer given by two one-electron orbitals with opposite spins and 
equal spatial wave function ip £ L^(M'^), ip being the unique minimizer of the restricted 
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Hartree functional <^ i— )• (/?), where £^ is defined in ([8]), below. In particular, 

if g = 2 and Z > 35 then mf£^^{ipi,ip2) > inf £^ {(p,ip). This statement is false for 
Z < 1.0268 (see the remark after Theorem II. 3p . which proves the necessity of some 
lower bound on Z other than Z > N — 1. 

The (unrestricted) A^-particle Hartree functional £^ is defined on arbitrary N tuples 
of functions ip^ G i?^(M'^) with J \ipi^{x)\'^dx = 1, and it is given by 



N 

E 

k=l 



+ 2 



2Z 



i<k 



Wi{x)?Wk{y)[' 

\x - y\ 



dxdy. 



(8) 



No orthogonality is assumed on (pi^ 



. . , (/9Ar. It is well-known that £^ has a minimizer 
if Z > — 1 and that minimizing orbitals are pointwise positive by a suitable choice 
of their phase [16]. To establish uniqueness of the minimizer we consider £^ as a 
spin-restricted Hartree-Fock functional with q = N, the spin-restriction being that 
(/)fc(x, s) = (pkix)6k,s- Then £^^{(f)i, ... , (f)^) = £^{^i, V'n), and the following 
theorem, with the exception of the bound on Z, follows from Theorem 11.21 with q = N. 



Theorem 1.3. Let N £ n and suppose that Z > 3^^(40A^ + 1GV2N - 8). Let 
(/?!,..., ifN be any minimizer of the Hartree functional ([8]). Then up to phases, ip^ = (p 
for all k where ip is the unique, positive minimizer of the restricted Hartree functional 
ip ^ £^{ip, ...,^). 
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Figure 1: The Hartree minimizer is unique for Z > Zc 



In general, the minimum of £ is strictly below the minimum of the restricted 
Hartree functional £^{ip,ip). In fact, Ruskai and Stillinger show that for Z < 1.0268 
and N = 2 the restricted Hartree functional is bounded below by —1, while for Z > 1 
it is straightforward to show that the unrestricted Hartree functional has its minimum 
below —1 [20j. A similar phenomenon of phase segregation is described in ^ for the 
two-electron Hartree functional with a confining external potential. 

We now briefly sketch our strategy for proving Theorem II. 11 First, we use the well- 
known fact that N functions Z^/'^ipi{Zx,a), i = 1, . . . , N , ipi £ H'^{R^ x {1, . . . ,q}), 
form a solution to the Hartree-Fock equations ([T]) if and only ipi, ... ,ipj\f is a solution 
to the new set of equations 

f . 2\ /^,2A/■ \ipk{y)\'^^i{x) - ipk{x)ipk{y)ipi{y) ^ 2^ . n 

V ~ R j '^'^^ z^J \^r^\ ^ " emix)- (9) 

These rescaled Hartree-Fock equations clearly exhibit the perturbative nature of the 
electron-electron repulsion in the large Z limit. Let Hp denote the rescaled Fock 
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operator on the left hand side in ([9]), see (fT3l) . the index P being the projection onto 
span{(/9i, . . . , 93jv}. For Z sufficiently large, Hp has exactly eigenvalues, counting 
multiplicities, near the eigenvalues —n^"^, /c = 1, . . . , s, of —A — 2/\x\. More precisely, 
if d is smaller than the gap nj'^ — {us + 1)^^ and if Z is sufficiently large, then the 
spectrum of Hp in Q := U|^-^[— n^^, — + d] consists of eigenvalues. The Hartree- 
Fock equations subject to (jH) are therefore equivalent to 

P = Xu{Hp). (10) 

Here we use that all shells are filled. By making Z even larger, if necessary, we can 
achieve that P i— )• Xn{Hp) becomes a contraction, and hence that (jlOp has a unique 
solution. This strategy obviously requires good eigenvalue estimates for Hp and good 
control of xn{Hp) — Xn{HQ) in terms of P — Q. Concerning the second point we 
develop a novel method that yields much better bounds than, e.g., a resolvent integral 
for Xn{Hp) would give. 

Our work was inspired by the paper of Huber and Siedentop on solutions of the 
Dirac-Fock equations |12j . Using the contraction principle, they solve an equation 
analog to (jlOp with the Laplacian replaced by the Dirac operator. We expect that 
our methods would allow to improve the results in |12] . Minimization problems for 
semi-bounded Dirac-Fock type functionals are studied in jlOj . and in the translation 
invariant case (no external potential) the minimizer is shown to be unique. When an 
external Coulomb potential is present the method of [lOj does not seem to work. For 
results on absence of a Hartree-Fock minimizer see [151 [21] . Existence and uniqueness of 
radial solutions to Hartree equations are derived in [19]. Uniqueness results concerning 
Hartree equations with attractive Coulomb interactions are established, e.g., in |13l [T], 
Several of the aforementioned results have been extended to pseudo-relativistic Hartree- 
Fock functionals and to Hartree-Fock functionals including a magnetic field [H El El [8] . 
Last but not least we should mention the fundamental paper of Bach on the accuracy 
of the Hartree-Fock approximation to the quantum mechanical ground state energy [3]. 

This paper is organized as follows. In Section 2 we introduce all notations and we 
prove eigenvalue bounds for Fock operators. Section 3 contains the proofs of all theo- 
rem in this introduction. In Section 4 we establish a theorem analog to Theorem 11.11 
in restricted HF-theory. Here, for proving existence of solutions we use the Schauder- 
Tychonoff theorem, which does not require as large values of Z as the contraction 
property does. 

Acknowledgement. M.G. thanks Alex Elgart, George Hagedorn and Timo Weidl for 
useful and inspiring discussions, Mathieu Lewin for explaining to him the uniqueness 
proof in [10], and I.M. Sigal for hospitality at the University of Toronto, where part of 
this work was done. F. H. is supported by the Studienstiftung des Deutschen Volkes. 
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2 Notations and Eigenvalue Estimates 



In this section we collect the operator- and eigenvalue estimates needed in later sections. 
We also introduce the definitions and notations used throughout the paper. 

The sets of bounded linear operators, of Hilbert-Schmidt operators, and of trace 
class operators in a separable Hilbert space % are denoted by C{H)^ 'T2{T-i) and 71 ('H), 
respectively. The corresponding norms are || • ||, || • II2, and || • ||i. Recall that 7i('H) C 
TiiV.) C CiH) and that ||K|| < ||K||2 < ||i^||i for all K G Ti{n). Let 

:= |P e Ti{L^{W^ X {!,..., g})) < P < 1, trP = iv}. 

The only reason for working with S^^q rather than with the set of self-adjoint projec- 
tions of rank A^, is that we need S^.q to be convex in Section HI 

For each P G SN,q there is a unique square-integrable kernel t = Tp such that 
Pip{x) = f t{x, y)ip{y)dy and there is a unique way to associate a density p = pp ^ 
with P, see Lemma ISTTl If P = Y2n^n\^n){v'n\, then 



T{x,y) = '^Xn^n{x)(Pn{y), (H) 
n>0 

p{x) = J2Xn\ipn{x)\'. (12) 

n>0 

For given P G Sjy^g we define a Fock operator Hp in L^(M'^ x {1, . . . , q}) by 

Hp = -A-V + ^{Up - Kp), (13) 

where V, Up are the multiplication operators associated with the real- valued functions 
V{x) = 2/\x\ and Up = p*V, and for ip € H'^, 

{Kp'ij:){x) := J V{x -y)T{x,y)'4'{y)dy. 

From the fact that V has a positive Fourier transform it is easy to see, using (jlip . that 
Kp > 0, and Up — Kp > by a straightforward computation. Hence Hp > —A — V> 
— 1. Moreover, by the Kato-Rellich theorem, Hp is self-adjoint on the Sobolev space 

//2(M3 X {1,...,^}). 

By remarks in the introduction we may consider as a function on SN,q(~^{P'^ = 
P} given by £^^{P) = £^^{ipi, ...,(/? at) where {ipi, . . . ,ipN) is any orthonormal basis 
of RanP. Explicitely, 



g^^{P)=tr[{-A-V)P]+^ J {p{x)p{y)-\P{x,yt)V{. 



X — y) dxdy. 



Lemma 2.1. Let P G SN,q- Then for all e > 

(i) V < -eA + e~^, 

(ii) Up < N {-eA + e-^), 
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1 1 / iV\^ 

(vii) -Up<e{-^-V) + -[e + -^ , 

Proof. Statement (i) follows from —A — V> —1 and from the scaling properties of —A 
and V with respect to the unitary transformation (/?(x, a) i— >• e~^^'^ip{'x/e, a). Statement 
(ii) follows from Up = p * V , from (i) and from J p[x)dx = N . 
To prove (iii), fix e > 0. By (i) and (ii), 

eV+^Up<[e + ^^ (-<5A + ri) 

for all > 0. Upon subtracting eV from both sides and making the choice 5 (e + ^) = e 
for 5^ the desired estimate follows. Inequality (iv) follows from — A < Hp-\-V and from 
(i). □ 

Proposition 2.2. Given iV G N, Z > and P G SN,q, let and denote the 
n-th eigenvalue, counting multiplicities, of the Schrddinger operator —A — V and the 
Hartree-Fock operator Hp, respectively. Then: 

(a) For all n^n, E^ <E^ and 

E^ <E^ + 2^ + 2^y^W+^- 

(b) If P minimizes the Hartree-Fock functional, then E^ obeys the following estimate: 

N-1 N-1 



E^<E^ + 2^^ + 2^^,/EfTl. 

Proof. From Up — Kp > 0, Kp > and Lemma l2.1l (iii) we see that 

1 _ If iV^ ^ 



-A - 1/ < < -A - 1/ + -Up < (1 + e) (-A - 1/) + - e + 

Z/ E \ /j 



for all £ > 0. By the min-max principle, this implies that 

E'^<El<{\^e)E'^ + -\e^ 



If N"-^ 



e\ Z 

Optimizing with respect to e yields the desired estimates of part (a). 

To prove (b) let . . . , Lpj^ be an orthonormal basis of RanP with Hp(pk = E^ (pk- 
Then the Hartree-Fock functional can be decomposed as 

S^^iVl, ■ ■ ■ ,<fN) = S^^i{ipi,...,(pN-l) + {<fN,Hp^N_iipN) 
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where and Hp^^-i denote the Hartree-Fock functional and Fock operator be- 

longing to the (A^ — l)-particle integral kernel TN-i{x,y) = 'Yl!k=i '^k{x)^k{y)- Since 
ijj I— )• £^^{ipi, . . . , ipN-i, ^) is minimized by ifN, 

E^= inf {i;,Hp^N-ii^)- 

■4>GD{Hp,N-i); 1^11=1 

ViGspan{i^i , . . . ,</3]V_ 1 }-'- 

Thus, using the min-max principle again, is bounded from above by the A^-th 
eigenvalue of Hp^jy^i. Part (a) applied to Hp^^i completes the proof of (b). □ 



3 Solving the Hartree-Fock Equations 

This section contains the proofs of all the theorems given in the introduction. For 
restricted Hartree-Fock theory see the next section. Our main tool is the following 
abstract result comparing the spectral projections xn{-^) and xn{B) of two self-adjoint 
operators A and B. 

Proposition 3.1. Let A, B : D C Ti ^ Ti be self-adjoint operators and let 0, be a 

bounded Borel set for which the spectra of A and B satisfy the gap conditions 

dist(a(A)nO,CT(5)\0) > 5, 

dist{a{B)r\n,a{A)\n)>5, ^^^^ 
for some 5 > 0. Suppose A and B have only point spectrum in Vt. Then 

\\xn{A) - xnm\2 < (iK^ - B)xnml + IK^ " B)xnm\iy^' ■ 

Proof. We only prove the proposition in the case where xniA) and xniB) are finite 
rank projections. The more general case of infinite point spectrum in $7 is left as an 
exercise for the reader. Let x(^) ■= Xn{A) and x{B) ■= Xn{B) for short. Then 
xiA)"^ = x(^) = x{^)* and similarly for x{B). It follows that 

IIX(^) - X{B)\\1 =tr(x(^)(l - xmxiA)) 

+ tTixiB)il-x{A))x{B)) (15) 

where we also used the cyclicity of the trace. To estimate the first term on the right hand 
side of ()15p we choose an orthonormal basis i<Pk)k=i of x(^) consisting of eigenfunctions 
of A: 

Aifk = Ek^Pk, £k G cr(^) n 17, k = l,...,n. 

By the gap assumption dHD, |A — efcl > 5 for A G (y{B) \ Q. and hence, by the spectral 
theorem, 

\\B - ek\>l - xn{B), k = l,...,n. 
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We conclude that 

n 

iT{x{A){l-xmx{A)) = Y.{vk,a-xiB))fk) 

k=l 

n 



fc=l 



k=l 



r^||(^-i?)x(^)||^. (16) 



The proposition follows from (I15p . (I16p and from an estimate similar to (I16p with A 
and B interchanged. □ 

Lemma 3.2. For all ip G H'^{R^ x {!,..., g}) and P,Q £ SN,q 

(i) \\{Up-UQM<4:\\P-QUVip\\ 

(ii) ||(i^p-KQ)<^|| <4||P-Q||2||Vv9||. 
Prooj. For all V G x {1, . . . , q)) 

\^{x)\ 



K^,(;7p-C/q)<^)| < 2 y" fiy|/5p(y)-/5Q(y)| ^ \^{x)\- 



dx 



\x - y\ 

< 4 1 dy\ppiy)-pQiy)\\\m'^^\\ 

< 4||P-Q||J|^||||V^|| 

by Cauchy-Schwarz, the uncertainty principle lemma \18\ section X.2] and by Lemma [5.11 
This proves (i). The proof of (ii) is similar. □ 

Proposition 3.3. Let P,Q £ SN,q H {P^ = P} and let C (— oo,0) be a bounded 
Borel set such that dist(Q, a{Hp) \ VL) > 5 and dist(il, (7{Hq) \ > 5 for some 6 > 0. 
Suppose moreover that tT{xn{Hp)) = N = tic{xn{HQ)). Then 



\xn{Hp) - xniHq) 
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< ^ (l + Vm) (\\^xn{Hp)\\l + \\V^xn{HQ)\\iy^' \\P - Q\\, 
<±(i + V2n)V2N\\P-Q\\,. 

The factor 8 in the last line may be replaced by 4 if both P and Q satisfy the Hartree- 
Fock equation (fTUj) . It may be replaced by 4-v/2.5 if P or Q satisfies ([TO]) . 

Remark. The set SN,qn{P^ = P} is a closed subset of the Hilbert-Schmidt operators 
on L^(M'^ X {1, . . . , g}) and hence it is complete with respect to the metric d{P, P) = 

\\P-Ph- 
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Proof. Applying Proposition 13.11 to Hp and Hq we see that we need to estimate 
\\{Hp — HQ)xn{Hp)\\^ and the same expression with P and Q interchanged. By defi- 
nition of Hp and Hq, 

Hp-Hq = Z-\Up - Uq) - Z-\Kp - Kq) 

and by Lemma [ 



\\{Up - UQ)xn{Hp)h < 4||P - Q\\i\\V^xniHp)\\^ 
\\{Kp - KQ)xn{Hp)\\2 < 4||P - Qh\\V^xn{Hp)\\^. 

The first inequality of the theorem now follows from ||P — < \/2iV||P — Q\\2- Here 
we use that P and Q have rank A'^. From Lemma l2.ll (iv) with e = 1/2 we see that 
-A < 2Hp + 4. Therefore 

\\V^Xn{Hp)\\l = tr {xn{Hp){-A)xn{Hp)) 
<AiTXn{Hp)=AN. 

If P solves the Hartree-Fock equation then ||\/— AxQ(-f^p)||2 = ^ A by the 

virial theorem and because —A — V> —1. □ 

For each R G S0(3) we define a unitary operator U{R) in L^(]R^) by 

[U{R)^]{x) ■.= ^{R-^x). (17) 

The following theorem describes the content of Theorem 11.11 in terms of rank N 
projections. 



Theorem 3.4. Under the assumptions of Theorem \l.l\ the Hartree-Fock equation (jlOp 
has a unique solution P £ Sj\f,g with the property that Z'^a{Hp\ P7i) is given by 
This P is of the form P = P' ® 1 with respect to L^(M^) where P' G Sjv/g^i and 

moreover P' = U{R)P'U{R)* for all R G 50(3). 

Proof. We first prove existence and uniqueness of P using Proposition 13.31 with 



fc = l 



-2 -2 , 



By Proposition 12.21 for Z large enough and all P G S^^q 

dist(0,a(//p) \^)> n-"^ - (n, + 1)-^ -^=:S, 

Z 

5 > 0, and trxn(-f^p) = A. Hence the Proposition 13.31 the remark thereafter, and the 
contraction principle imply that the equation P = xn{Hp) has a unique solution P 
provided that 

-I (l + V2n) V2N < 6. (18) 



Z 

This is satisfied for Z obeying 
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We next show that P is of the form P' 1 in L'^{R^) ® C where P' G 5^^/,^!. To 
this end we consider the modified Hartree-Fock equation 

P' = Xq(4'') (19) 

where 

H^^) ■.= -A-V+ ^Up, - ]-Kp, 

in L^(M'^). The eigenvalues of iJp,^ satisfy the estimate given by Proposition 12.21 a). 
Therefore, the arguments above show that (fT9]) has a unique solution P' because ([T8|) 
holds by assumption ([3]). From Hp'(^i = Hp) 1 and (fT9]l it follows that P' (g) 1 solves 
the Hartree-Fock equation (fTOjl . Hence, P = P' 1 by the uniqueness of the solution 
to PU]). 

Finally we prove that P' commutes with U{R). From the spherical symmetry of V 
it follows that U{R)H^^}U{R)* = where P'{R) = U{R)P'U{R)*. Using (fT9]l. 

we conclude that 

P'{R) = U{R)xn{H^^hu{Rr = xn{y{R)H^^}u{RY) = 

which implies P'{R) = P' because the solution to (|19p is unique. □ 

Proof of Theorem Theorem 13.41 tells us that P = P' ® \ where P' commutes 
with all rotations. Moreover, Hp = H^p} (^1, P' = xn{Hp)) and Hp) commutes with 
all rotations as well. Let S be an eigenvalue of Hp) . The eigenspace associated 
with Efc carries a representation of SO{3) given by p!7|) . Its irreducible subspaces are 
spanned by functions of the form 

|xrV(|x|)l;r„(x), / e L2(M+), (20) 

where Yi^ denotes a spherical harmonic. Now fix / and m and let Him denote the 
space of all functions of the form (|2Up with arbitrary /. This space is reducing for 
both —A — V and Hp) . Since the spectrum of (—A — V) \ Him in (-co, 0) consists of 
the simple eigenvalues — 1/n^, n > Z + 1, it follows from Proposition 12.21 and from the 
assumption on Z that Hp) \ Him has exactly one eigenvalue in each of the intervals 
[— n^^, —{uk + 1)~^) with nfc > / + 1. This completes the proof of Theorem II. li □ 

Proof of Theorem For Z > N — 1 the Hartree-Fock functional is known to have 
a minimizer and any minimizer P is the spectral projection onto the spectral subspace 
of Hp associated with the lowest N eigenvalues [16]. li N = q Yln=i ""^ some s, then, 
by Proposition 12.21 the lowest N eigenvalues of Hp belong to [—1, — (s + 1)~^) provided 
that 4(A^ — 1)/Z < — (s + 1)~^. Thus for sufficiently large Z Theorem 11.11 implies 
uniqueness of the minimizer of the Hartree-Fock functional as well as the assertions on 
the one-particle orbitals. The bound on Z is obtained by inspection of the proof of 
Theorem 11.11 keeping in mind that P is a minimizer which is given. Hence the improved 
bounds from Proposition 12.21 and Proposition 13.31 are available. □ 
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Proof of Theorem[TIE Given 931, . . . , c^at E with / \ipk{x)\'^ dx = I let $1, . . . , ^>Ar e 

H^{R^ X {1, ... , N}) be defined by 

$fc(x,s) := (/9fc(x)4^, k = l,...,N. (21) 

Then <i>i, . . . , <I>Ar are orthonormal in L^(M^ x {1, . . . , A^}) and 

S^^i^i, ...,'^n) = S^i^u ■■■,^n) (22) 

by the definitions of and £^ . By (j22p . the minimization problem for £^ is 

equivalent to the minimization problem for with q = N in the restricted class of 
orbitals of the form (121 p . By Theorem ll.2| the (unrestricted) Hartree-Fock functional 
for q = N and Z sufficiently large has a unique minimizer P £ Sn,n, which is of the 
form P = P' 1 with respect to L'^{R^ x {1, . . . , N}) = L'^{R^) C^. P' here is 
a rank one projection that commutes with rotations. Hence (P'lp) = ^{f, ip) with a 
spherically symmetric function ip. It follows that £^^{P) = £^ {(p, . . . ,(p). It remains 
to determine the condition on Z for uniqueness of the Hartree-Fock minimizer in the 
present case where q = N. To this end we use the improved eigenvalue estimate from 
Proposition 12. 21 b). The gap condition becomes 

^■=1- > (23) 

and the contraction condition reads 

4 (1 + Vm) V2N < S. (24) 



z 

Both, ([23]) and ([MD are satisfied if 



z>i2iv + 15v2iv-5. 

3 3 3 



□ 



Corollary 3.5. Under the assumptions of Theorem \3.4\ let Pz be the unique solution 
provided by this theorem. Then 

JimP. = tx{_„,.}(-A-y). (25) 
fc=i 

Proof. Let Pqo denote the right hand side of ()25p . A copy of the proof of Proposition [33] 
with P = Pz and Hq replaced by i/oo = —A — V shows that 

\\PZ - Pooh ^JzV + ^) ^11^^112 = V ^ ^) ^ ° ^ 

Here, by Proposition l2.2| 5 can be chosen independently of Z for Z sufficiently large. □ 
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4 Restricted Hartree-Fock Theory 

In this section the Hartree-Fock functional is restricted to one-particle orbitals of the 
special form ([5]). This will allow for all the electron configurations found in noble gas 
atoms. In this section we set g = 1 to simplify the presentation. 

By a shell index we mean a pair of integers (n, tj with n > 1 and < £ < n — 1. 

Theorem 4.1. Let {ni,ii), . . . ,{ns,is) be given, pairwise distinct shell indices with 
ni < n2 < ■ ■ ■ < Ug. Suppose N = X]j=i(2^i + 1) ^''^^ let := n~'^ — {ug + 1)~^. 

(i) If Z > an/As then there exist normalized functions fi, . . . , fg G L^(IR4-) such 
that the N functions 

'Pjm{x) ■■= -^^fj{\x\)Yi^jn(.x), j = l,...,S, TU = -Ej . . .£j, 

solve the Hartree-Fock equations ^ with eigenvalues Sj satisfying 

1 1 4iV 

-^<e,<-^ + —. (26) 
n. n- Z 

(a) If Z satisfies ([3]) then the functions fj in (i) are unique up to global phases. 

Remark: A result similar to part (i) of this theorem with the weaker assumption 
Z > is described in Section III. 3 of |17) . However, Lions' argument is based on the 
unproven assertion that all eigenvalues of a radial Hartree-Fock operator are simple. 

To prove Theorem 14.11 we solve the fixed point equation 

s 

P = F{P):=Y,XnAHp)7ri^ (27) 
i=i 

where tti denotes the orthogonal projection associated with the eigenvalue i{i + 1) of 
the square of the total angular momentum operator, and 



1 1 AN 

nf n] Z 



The spherical symmetry will be imposed by restricting F to the subset 
S'^"^ := {P e SN,i\U{R)PU{Ry = P, for aU R € 50(3)} 



of Sn,i- To find a solution of (j27|) we use the Schauder-Tychonoff theorem. Its unique- 
ness will follow from the contraction principle. 



Lemma 4.2. Suppose the hypotheses of Theorem \4- 1\ are satisfied and that 6 := Ag — 
AN/Z > 0. Then FiS'j^"") C S'^"" and for all P,Q e S'^"" , 

||F(P)-F(Q)||i<32^||P-Q||i. (28) 



lfP,Q£ S'y"'n{P^ = P} th 



en 



\\F{P) - F{Q)\\2 < J^V2N{1 + V2N)\\P - Qh- (29) 
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Remark. The Lipshitz constant in (I29p agrees with the one in Proposition 13.31 

Proof. We first show that F(S'^™') C S"^™. For P G -S^™, the Fock operator Hp 
commutes with all rotations U{R) and hence so does its spectral projection XQj{Hp). 
It follows that Fj{P) := xnj{Hp)7ri. = TTi.XfijiHp) is an orthogonal projection and 
that 

FjiP)FkiP) = 6jkFjiP) (30) 

because {nj,£j) ^ (jikAk) for j ^ k. Hence F{P) = Ylj=i ^j(^) orthogonal 
projection that commutes with all rotations. As for the trace of F{P), we note that 
Hp \ TTi-'H has exactly 2ij + 1 eigenvalues in ilj, counted with multiplicities. This 
shows that trF(P) = Ej=i tr -Fi(-P) = Ej=i(2^i + 1) = A^- 

By inspection of the proof of Proposition 13. H using (f30]l we see that 

\\F{P) - Fm\l = tr F(P)(1 - F(Q))F(P) + trF(Q)(l - F(P))F(Q) (31) 

where 

s 

trF(P)(l - F{Q))F{P) = ^trF,(P)(l - F(Q))P,(P) 

s 

<r2^||(i7p-/7Q)P,-(P)||i. (32) 
i=i 

Here it is important that Hp and Hq are considered as operators on Tri.L'^(W^), which 
have 2ij + 1 eigenvalues in fij separated by a gap of size 6 from the rest of the spectrum. 
To bound ||(^^p — HQ)Fj{P)\\2 we use estimates from the proof of Proposition 13.31 as 
well as IIP - QII2 < IIP - Qlli. We find 

\\{Hp - HQ)FjiP)\\2 < ^y2Z7TT||P - Qlli. (33) 

The proof of ([28]) is now completed by combining ([3T]) , ([32]) , ([33]) with 

\\F{P) - P(Q)||i < Vm\\F{P) - P(Q)||2 (34) 

and with estimates similar to ([32p . ([33p where the roles of P and Q are interchanged. 
(I34p follows from the fact that P(P) and F{Q) have rank N. For proving (I29p we use 
||P-Q||i < \/2iV||P-Q||2 instead of ([MD- The rest is similar to the proof of (l28]) . □ 

Proof of Theorem \4-l\ To prove (i) we check that the map F : — )• S*^™" sat- 
isfies the hypotheses of the Schauder-Tychonoff theorem. By construction, S^™" is 
convex and closed in 71(L^(M^)), and F : — )• is continuous by Lemma H^l 

To prove compactness of the map F let (P„) be any sequence in S^™. Then F{Pn) 
may be written in the form 

N 

-f'l-^nj - 2^l93fc I, ^Pn93fc -efc 

k=l 
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with eigenvalues e^"^ < — (n<j + 1) ^ for all k,n. By Lemma \2A] the sequences )n 
are bounded in H^{R^). We may assume they are weakly convergent in and hence 
locally convergent after passing to a subsequence of (Pn)- Hence F{Pn) will be con- 
vergent provided that ip^j^\x) — ?• as |x| — t- oo uniformly in n. To prove this uniform 
decay we pick a function x ^ C°°(M; [0, 1]) with = for i < 1 and x(0 = 1 
t > 2. Let xr{x) = xi\x\/R) for R>0. The IMS-formula tells us that 

2xRiHp„ - el"^)xR = xUHp„ - 4"^) + {Hp„ - e^^^)x% - [[Hp„,XR],XRi 

where [[Hp^,xr],Xr] = 2|Vxi?P - Z-'^[[Kp^,xr\,Xr] and 

Xr{Hp^ - et^R > {-2R-' - et^)xl > cxl 



for some c > if i? > iio and Rq is large enough. It follows that 

' U^\[xrAxr,KpJ]^^-^) +0{R^ 



cWXR^l 'II' - 2Z 



= 0{R-'), (R^oo), 

uniformly in n. In the last step we expanded the double commutator, we used | r^") {x,y)\ < 
p(n)^2;)i/2p(")^y)i/2^ Cauchy-Schwarz, the spherical symmetry of p^"^ and Newton's the- 
orem to replace \x — y\~^ by min(|x|~^, Here r^") and /o^"^ denote the kernel 
and density of P„. Statement (ii) in Theorem 14.11 follows from (|29p and the remark 
following Lemma |4.2[ □ 

5 Appendix 

The following lemma seems to be well-known, but we have not found it in the literature. 
Lemma 5.1. There exists a bounded linear mapping 

p:ri(L2(M")) ^ L\W) 

K ^ PK 

which is uniquely determined by 



\'lp){(p\ H> 'll^{x)(f{x). 
Furthermore, p has the following properties: 

(i) J Pk{x) dx = iv{K) 

(ii) j \pK{x)\dx <t'£\K\ 

Proof. The map p is unique because, by linearity and continuity, it is completely deter- 
mined by its action on operators of rank one. For the proof of existence let J G C^(M") 
be a real, non-negative and even function with J J{x) dx = 1. Let Je{y) '■= e~"J(y/e) 
and Je,xiy) ■= Jeiy - x). Then 

{Je,x,^) = / Js{x-y)ip{y)dy = {Js*(p)ix). (35) 
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For given K E Ti{L\W)) we define pk G L^{W) by 

PK ■■= - lim pK,e 
e—^0 

PK,e{x) ■■= {J£^x,KJe^x)- 

We claim that the limit ()36p exists and that 

n>0 

if K = Xn\ipn) {lpn\ 



(36) 



(37) 



n>0 



denotes the singular value decomposition of K, that is > and (pn, ipn are orthonor- 
mal bases. Using ([35]) . (j371) . ||Je * (/?|| < || j£||i||c/?|| = \\ip\\ and X]„ A„ = tr|i^|, we see 
that 

ll/OK.elli < tr|Er|, forane>0. (38) 

Hence if the limit (j36p exists, then (ii) holds and K i— )• /jj^ is a continuous linear map. 
Moreover, if denotes the A^-th partial sum of ([37|) . then, by ([38|) . 



< tr - Kjvl ^ 0, (iV ^ oo), 

uniformly in e > 0. Therefore it suffices to prove the existence of pxf^^ which follows if 
we prove that p maps the rank one operator to p{x)ip{x). Indeed, 



(Je * 'p){J£ *tlj) - iflp 



{Je*p- p)lp + {Je * p) \ {Je * ^) ~ „^ 
< II Je * — IIV'II + II * l/^ll II * V' — V'll — ^ (ff — )• 0) 

where we used that Je *(/?—>• in L-^(R"). It now remains to prove (i). This follows 
from 



/ pK{x)dx = y2i^n / Pn{x)'4)n{x)dx = y2^n{'^n, 
n>0 •' n>0 



'Pr. 



Pn 



n>0 



n>0 



□ 
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